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The effect of boundary deformation on the non-separable entanglement which appears in the classical elec¬ 
tromagnetic field is considered. A quantum chaotic billiard geometry is used to explore the influence of a 
mechanical modification of the optical fiber cross-sectional geometry on the production of non-separable entan¬ 
glement within classical fields. For the experimental realization of our idea, we propose an optical fiber with a 
cross section that belongs to the family of Robnik chaotic billiards. Our results show that a modification of the 
fiber geometry from a regular to a chaotic regime can enhance the transverse mode entanglement. Our proposal 
can be realized in a very simple experimental set-up which consists of a specially designed optical fiber where 
non-entangled light enters at the input end and entangled light propagates out at the output end after interacting 
with a fiber boundary that is known to generate chaotic behavior. 


I. INTRODUCTION 


Recent research in quantum chaos has uncovered the rela¬ 
tionship between chaos and entanglement both theoretically 
and experimentally M. It is well known that classically 
chaotic trajectories exhibit rich topological structures in the 
classical phase space jib- Recently, there has been an in¬ 
creasing interest in relating quantum entanglement and clas¬ 
sical geometry. In this article, we explore how quantum¬ 
like entanglement can be intimately connected to the geom¬ 
etry of the domain boundary within which the wavefunc- 
tions are confined. In this context, a quantum chaotic bil¬ 
liard in the non-relativistic quantum regime is best suited 
for this study, with the effects of spatial geometry on entan¬ 
glement and its connection to Hamiltonian chaos to be an 
as yet unexplored direction. Classical chaos has been ob¬ 
served in very simple systems like the Robnik billiard |0] 
and the Bunimovich stadium|01. It is already known that a 
quantum system in a classically chaotic regime can generate 
higher entanglementlfH [sj- fl^ . In terms of its implementation, 
a recent experiment has demonstrated the presence of quan¬ 
tum chaos in light beams m. While quantum entanglement 
has been practically implemented in sophisticated discrete- 
variable systems such as an NMR based system ifT^ . there is a 
recent interest in the application of continuous variable entan¬ 
glement for the purpose of quantum computation fil . The 
approach based on the analogies between quantum physics 
and optics via the classical electro mag netic waves is partic¬ 
ularly relevant in this context as in |[l^. It was Spreeuw who 
realized the existence of non-quantum optical entanglement 
sm. According to Ghose, the EPR like correlation can 
also appear in classical electrodynamics |23t]. More re¬ 
cently it is shown that the classical electrodynamics contains 
a Hilbert space structure which makes it similar to quantum 
mechanics in certain aspects ll23l . In the recent experiments, 
researchers have explored the non-quantum entanglement in 
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the classical optical field 1^ . It is well known that the 
Maxwell’s equations of the transverse modes of the electric 
field in an optical fiber can be transformed into an Optical- 
Schrodinger equation, where the z-axis of the optical fiber is 
analogous to time ll27l - [^ . Hence, the classical electromag¬ 
netic wave equation can be used to mimic the quantum fea¬ 
tures and facilitate a proper experimental implementation. H. 
J. Stockmann et al. has already explored this type of anal¬ 
ogy between electromagnetic waves and quantum mechanics 
to explore qua ntum chaos using a two-dimensional microwave 
cavity Here, we use geometries inspired from stud¬ 

ies of quantum chaos to explore the geometric dependence of 
entanglement. We introduce the idea of exploiting an opti¬ 
cal fiber with a core cross-section that has the geometry of 
a chaotic billiard to create a highly entangled light beam. 
More specifically, we analyze in this paper an optical fiber 
which possesses the cross sectional geometry of the Robnik 
billiard. By analyzing the propagation of a classical electro¬ 
magnetic wave in these fibers, we can explore the geomet¬ 
ric dependence of the non-separable entanglement in terms of 
the lowest-order eigenmodes, as well as using coherent and 
squeezed coherent states as the initial wavepackets. 


II. TWO-DIMENSIONAL OPTICAL-SCHRODINGER 
EQUATION 


It is well known that a light beam propagating in an optical 
fiber can be described by the following equation. 
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which is a two-dimensional Optical-Schrodinger equation 
with an optical-wavefunction in the transverse x and y vari¬ 
ables. The longitudinal coordinate z along the fiber axis is 
similar to time in the standard Schrodinger equation & 
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[Ij]]. This type of analogy between quantum mechanics 
and wave optics has been also pointed out in many other works 
lIsB-S- Here no(z) is the refractive index along the beam 
axis and is the wavelength of the electromagnetic wave. 
Hence, the optical quantum Hamiltonian can be written as, 


H = 


1 

no(z) 



U(x,y,z) 


( 2 ) 


and the optical potential function has the form U(x, y, z) = 
2 nl{z) [^o^(-g) “ n^(x, y, z)]. This type of two-dimensional 
Hamiltonian has already been studied in the context of chaotic 
ray propagation in an optical fiber . In the usual sense, the 
classical limit corresponds to ?i —0, while in the Optical- 
Schrodinger equation this limit corresponds to Au, —?> 0 which 
is just the classical ray optics approximation. By properly 
choosing the refractive index profile, we can define our po¬ 
tential function. 

Let D be the domain of the fiber core cross-section in the 
two-dimensional xy plane. The optical fiber is a dielectric 
medium with refractive index n{x,y,z). It is known that for 
a step-index fiber, n{x, y, z) = uq in the core with (x, y) G 
D U dD and n{x, y, z) = ni in the cladding if (x,y) G D. 
We assume that the outer region of the cladding is sufficiently 
large so that it can be taken to be infinite when considering 
guided light in the core and near the core-cladding boundary. 
Under this situation the components of the monochromatic 
electric and magnetic fields obey the Helmholtz equation. In 
order to study the time evolution of the wavepacket, we define 
our two-dimensional potential as a function which is zero in¬ 
side the domain (fiber core region) and relatively high outside 
the domain (cladding region). Hence our problem can be rep¬ 
resented by the Helmholtz equation in a 2D domain which has 
been employed to describe diverse quantum chaotic billiards. 
In other words, our problem of the wavepacket propagation in 
an optical fiber has been reduced to a quantum chaotic billiard 
problem. The two-dimensional potential is taken as zero in¬ 
side the domain and a high value of around 10^^ outside the 
domain. This is enough to make the wavefunction to vanish on 
the outer boundary. In order to represent the two-dimensional 
potential on a grid, the Robnik billiard is represented in carte¬ 
sian coordinates which is given by 

((m -f Bf +v^- 2B{u + B)f = A^{{u + Bf + v^). 

(3) 


Using the cartesian form, the region inside and outside the 
Robnik billiard can be easily determined. Once the two- 
dimensional potential is obtained numerically, the wavefunc¬ 
tion is propagated via the split operator method. Then, the 
density matrix is computed to evaluate the von Neumann en¬ 
tropy of entanglement Sm- 

In the classical limit, we study the dynamics of the classical 
billiard via ray optics while the analogous quantum mechan¬ 
ics represents the wave optics description in terms of the elec¬ 
tromagnetic wave and its propagation in the waveguide. The 
quadrature components x, y, Px and py are the position and 
momentum coordinate in the plane perpendicular to the fiber 
axis. This is clarified in Fig.[T]where we have shown the wave 



Figure 1. This figure shows the geometry at the core of a family of the 
Rohnik optical fiber. When A = 0, the core of the optical fiber takes 
the standard circular shape. By increasing the deformation parameter 
A, different quantum chaotic optical fibers are obtained. A schematic 
of the optical ray propagating in the ^-direction is also shown as well 
as its projection on the xy plane which contains patterns reminiscent 
of a quantum chaotic billiard. Note that the 2 axis is the analog of 
time in the standard Schrodinger equation. 


propagation in the optical fiber and its projection on the (x, y) 
plane. It is observed that the projection on the (x, y) plane 
is a quantum chaotic Robnik billiard. Note that in our case 
X, y, Px and py are exactly the coordinates and momentum 
variables of the quantum chaotic Robnik billiard. In the ray 
optics limit, the path of light rays illustrate the chaotic proper¬ 
ties while in the wave optics regime, the interference of light 
waves gives rise to quantum-like effects. 


III. MECHANICAL MODIFICATION OF THE QUANTUM 
ENTANGLEMENT 

The goal of this article is to elucidate the effect of the do¬ 
main boundary geometry on non-separable entanglement. Our 
results would illustrate the possibility of manipulating optical 
entanglement via modifying the mechanical design of the core 
geometry. We shall entirely focus on the Robnik billiard due 
to its relative simplicity in terms of numerical computations. 
Our idea can be intuitively grasped from the following fact 
that the eigenfunctions of a rectangular billiard is formed by 
a tensor product of x and y modes, and hence it is not en¬ 
tangled. By slightly deforming the rectangular geometry to 
a Bunimovich stadium, the eigenfunction would be changed 
entirely. The new eigenfunctions of the Bunimovich stadium 
then becomes a linear combination of the eigenfunctions of 
the rectangular billiard with proper boundary conditions and 
it is no longer a product state. The same is true for the Robnik 
billiard which results from deforming a circle. The eigenfunc¬ 
tions of the circular billiard 10 can be expressed as a tensor 
product state in terms of the r and 0 variables as follows: 

(4) 

In the radial direction, the eigenmodes are described by the 
Bessel functions Jmijmnx), while the polar coordinates give 






3 


the wavefunction 0 = 6+*"*® and jmn describes the nth zero 
of the Bessel functions of order m. Hence, there is no en¬ 
tanglement between the polar and the radial part of the wave- 
function for the circular geometry. By deforming the geome¬ 
try, quantum entanglement starts to appear. For example, the 
eigenfunction of the Robnik billiard can be expanded in terms 
of the eigenfunctions of the circular billiard as follows; 

N 

i,{v,e) = ^ (5) 

m—0 

The value of /3 can be determined by minimizing the norm 
of the wavefunction on the boundary, which then gives the 
complete information about the wavefunction ip{r,6). This 
technique is often used in the numerical implementation of 
the method of particular solutions EHH. Equation |5] shows 
that the wavefunction 9) is in the form of the Schmidt 
decomposition, where the Schmidt basis are the eigenfunc¬ 
tions of the circular billiard and the entanglement can be de¬ 
termined easily from the coefficients Cm- As the deforma¬ 
tion increases, we need more terms to approximate the eigen¬ 
function. Hence, the entanglement entropy would increase. 
When the beam paths configuration becomes highly chaotic, 
a larger number of Schmidt modes are required to approxi¬ 
mate the eigenfunction which leads to a higher entanglement. 
This results from the fact that the field modes start to interact 
through the special geometry of the deformed circular bound¬ 
ary. This simple observation has far reaching consequences 
since it implies the possibility of manipulating entanglement 
via mechanical methods which is technically easy to achieve. 
In addition, quantum chaos would serve to indicate the critical 
geometry that gives a maximum entanglement in the system. 

In order to study the continuous variable entanglement of 
the transverse mode '0(r, 9) in quantum chaotic optical hbers, 
we hrst obtain the reduced density matrix of the hrst subsys¬ 
tem Pi by integrating over the field variable 9 in terms of the 
bipartite wavefunction 9) as follows: 



(a) 




71 = 1 



(b) 


Figure 2. The von Neumann entropy of entanglement of the four 
lowest eigenfunctions vs the deformation parameter A of the Robnik 
fiber is shown in (a), while (b) shows the four lowest state eigenfunc¬ 
tions with A = 0.499. It can be seen that the von Neumann entropy 
of entanglement Svn in the eigenmodes saturates as the geometries 
of these fibers approach the completely chaotic regime. The solid 
black curve in (a) shows the average von Neumann entropy of en¬ 
tanglement of the ten lowest eigenmodes Svn for different boundary 
geometries of the Robnik fiber. 



pi(r,r)= J tj;{r,9)ilj*{r ,9)VT^d9 . (6) 

To quantify the continuous variable entanglement, we com¬ 
pute the von Neumann entropy of the entanglement by using 
the numerical methods proposed by Parker et al. 11^ and 
Bogdanov et al. a. The von Neumann entanglement en¬ 
tropy of the reduced density matrix is given by, 

Svn (0 = - X! (7) 

where pi are the eigenvalues of the Hermitian kernel pi (r, r ) 
as shown in Eq. [^below. The eigenvalues are computed from 
the Eredholm type II integral equation of pi(r, r ), which is 
given by 

y pi(r,r')(/)i(r')dr' = rji ^^(r), (8) 

with (pi (r) being the corresponding Schmidt eigenfunction. 


IV. ROBNIK FIBER AND ENTANGLEMENT 

The cross section of a Robnik fiber has the shape of the 
Robnik billiard which is shown in Eig. [1] As light passes 
though the Robnik fiber the transverse modes of the fiber get 
entangled. We are interested in the optimum geometry which 
gives a high entanglement in the system. The cross sectional 
geometry of the Robnik fiber is a conformal transformation 
w = AZ{9) -F BZ{9)^ of the unit circle Z{9) on the complex 
plane. According to Robnik |@], the cross-sectional geometry 
can be written in the parametric form as follows, 

u = A cos 9 + B cos 29 (9) 

u = Asin0-F i3sin20. (10) 

The area of the billiard is given by Area = 7r(Fl^ -F 2B^) and 
by means of a reparametrisation, A and B can be written as 
follows: A — cosp and B — sinp, where p is given by 

p = tan“^ (Av/2)- This allows us to define a family of closed 
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(a) 



Figure 3. A plot of the entanglement of the squeezed coherent state 
versus the deformation parameter A of the Rohnik fiber is shown in 
(a). We observe that the entanglement is higher when the system 
is in the chaotic regime and the initial squeezing has enhanced the 
entanglement in the system. In (b), we have plotted the entangle¬ 
ment dynamics as the light propagates along the 2 :-axis. It is impor¬ 
tant to note that the light enters the fiber without any entanglement. 
Entanglement increases as it propagates, and as it traverses a dis¬ 
tance of 50 units the entanglement saturates. Note that the blue and 
green curves represent the coherent state and the squeezed coher¬ 
ent state respectively. We have chosen the initial state at the point 
{x,Px,y,Py) = (0.25,0.1,0.0,0.1). 


analytical boundaries of constant area, Area = tt, which de¬ 
pends on the values of the deformation parameter A. Rohnik 
had shown that the parameter A can be chosen within the in¬ 
terval 0 < p < Psing = tan“^ ("Tf ^ smaller values of A, 
the billiard exhibits a regular dynamics since it is closer to the 
circular geometry. An increase in the parameter A shifts the 
system from regular to chaotic. We analyze the dependence 
of the quantum-like entanglement as we vary the deformation 
parameter A from the regular to the chaotic regime. 


A. Entangled Rohnik eigenmodes 

We have varied the fiber cross section by adjusting the de¬ 
formation parameter A. As the boundary geometry changes. 


we compute the lowest order eigenfunctions, so that the en¬ 
tanglement in the radial r and the angular 9 modes can be 
evaluated. Based on our earlier discussion, the entanglement 
in the circular geometry between the r and 9 variable is zero 
since the eigenfunction is a product state. As we change A, the 
quantum-like entanglement in the system starts to increase. In 
Fig. Ilia), we have plotted the entanglement of the eigenmodes 
against different values of the deformation parameter A. In the 
same figure, we have also displayed the corresponding bound¬ 
ary geometry. We observe from Fig. ^a) that the entangle¬ 
ment dynamics increases as A increases. As the boundary ge¬ 
ometry leads to chaotic behavior, the entanglement in eigen¬ 
modes increases even further and as the system approaches 
the highly chaotic regime the entanglement entropy saturates. 
Thus, the degree of chaos in the system is directly linked to 
the level of entanglement in the light pulse albeit there is a 
maximum attainable level in lieu of the eventual saturation. 
We have plotted the eigenmodes corresponding to the satu¬ 
ration level at A = 0.499 in Fig. Ilfb). From Fig. IS a), we 
observe that the lowest energy mode n = 0 has higher en¬ 
tanglement compared to the excited mode n = 1. This can 
be understood by examining Fig.|2jb), where the ground state 
eigenmode is seen to be more severely affected by the defor¬ 
mation of the boundary, while the first excited mode n = 1 
contains a nodal line through its center and has a form that is 
closer to the eigenmodes of the circular cross-section. Hence, 
the deformation in the eigenfunction is smaller for the n = 1 
mode. The number of basis terms to describe this eigenfunc¬ 
tion is smaller, which implies that its entanglement is smaller. 
The same argument applies to the excited modes n = 2 and 
n = 3, respectively. These results indicate that the eigen¬ 
functions that are affected more by the boundary deformation 
would have a higher entanglement. It thus provides insights 
into how quantum chaotic geometry can be effectively used to 
produce highly entangled states. 

B. Effect of squeezing and entanglement in the Rohnik fiher 

While it is typical to propagate a Gaussian beam in an opti¬ 
cal fiber, one may wonder the effect of initially squeezing such 
a beam before launching it into the fiber. This question has led 
us to investigate the difference in entanglement dynamics with 
respect to the initial coherent and the squeezed coherent state 
as the light wave propagates in the horizontal z direction. For 
our investigation, we use the Hollenhorst and Caves definition 
of the squeezing operator 1431 - 1?^ and we employ the coor¬ 
dinate representation of the squeezed coherent state ll4^ for 
the numerical computations. The squeezed coherent state is 
defined as 

|afe,a) = C(afe)5(Cfc)|0), (11) 

where the displacement operator D and the squeezing opera¬ 
tor S is given by 

D{ak) = exp (afettfe''' - alcik) , (12) 

5(Cfe) = exp i^Ckdk'^ - ^Ck*dk‘^) , (13) 
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Figure 4. The probability densities of the entangled coherent states coming out of the Robnik fiber with length z — 50 units 

are shown in (a), (b) and (c) with A = 0.15, A = 0.2 and A = 0.499 respectively. The probability densities are shown in the 

cross-sectional u — v plane of the fiber. The blue color shows the regions where the probability is zero and the red color shows the regions 
of higher probability. The entangled state is obtained from a tensor product coherent state centered at {x,px,y,py) = (0.25, 0.1, 0.0, 0.1) 
with Au, = 0.01. The entangled state occurs after the light beam has propagated a distance of z = 50 units. Note that we have also indicated 
the von Neumann entanglement entropy of the polar coordinates below these figures. In (d), (e) and (f) the corresponding classical phase 
space is shown for A = 0.15,A = 0.2 and A = 0.499 respectively. From the figure it can be observed that the time evolved coherent state has 
a higher entanglement when the corresponding geometry possesses a larger degree of chaotic behavior. 


and ak = (^k = \fk\ are complex numbers 

and ak are related to the phase space variables {qk,Pk) in the 
following manner 

ttfc =+ ipfc), (14) 

with k = 1,2, respectively. According to Mpller et al. 

the squeezed coherent state in the position basis can be written 

as 

V'(x,afc,C/c) = (;^)^^‘^(coshrfe+ e‘®sinhrfe) 

exp {-^ (X - qif + ^Pi(x - gi/2)} 

(15) 


The tensor product state of this wavefunction is used to study 
the quantum-like entanglement dynamics for different squeez¬ 
ing parameter values. 

In Fig. [3ja) the entanglement maximum Sm is plotted 
against different fiber geometries. The entanglement maxi¬ 
mum Sm is taken from the von Neumann entropy of entan¬ 
glement Svn{z) which is computed for the propagation of 
the initial states along the z-axis with a fixed geometry. The 
initial tensor product coherent state and the squeezed coher¬ 
ent state is propagated for a different boundary geometry and 
we take these initial states to be centered at {x,Px,y,Py) = 
(0.25, 0.1, 0.0, 0.1). We have computed the quantum-like en¬ 
tanglement present in the radial and polar variables of the op¬ 
tical wavefunction. In Fig. [3b) the entanglement dynamics 
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of the coherent state and squeezed coherent state is plotted 
against the z-axis. When the optical wave with a tensor prod¬ 
uct coherent state enters the fiber (at z = 0), its entanglement 
is zero. During its propagation along the z-axis the trans¬ 
verse modes get entangled and after traveling a distance of 50 
units along the z-axis, the entanglement becomes saturated. 
In other words, we observe that a non-entangled light enters at 
one end of the fiber and comes out entangled at the other end. 
In the case of a highly quantum chaotic fiber with A = 0.5, this 
entanglement can yield a higher value. Note that we have cho¬ 
sen the initial condition {x,Px,y,Py) = (0.25,0.1,0.0,0.1) 
because it lies entirely inside the domain for the chosen pa¬ 
rameter range of A. The wavelength (analog of K) given in 
Eq.[T]is selected as A^, = 0.01. The value of A^, is set small 
enough in order for the wavepacket to be contained entirely in¬ 
side the defined two-dimensional cross-sectional domain. We 
have also used the equally position squeezed Ci = C2 = C 
wavepacket at a relatively small value of ^ = —0.25 due to 
the numerical considerations. From Fig.[3we can clearly see 
that as the core boundary deformation A increases, the produc¬ 
tion of entanglement is found to increase for both the case of 
initial coherent state and initial squeezed coherent state. This 
is again observed in Figs.|4ta),|4fb) and|4jc) according to the 
corresponding von Neumann entropy of entanglement. It is 
well known that when the deformation parameter approaches 
A = 0.5, the dynamical behavior of the classical system is 
highly chaotic (see Figs. Ud), |4je) and|4jf)) with the higher 
entanglement maxima observed in the corresponding quan¬ 
tum system. Thus, our results clearly indicate the effect of the 
geometrical deformation on entanglement dynamics of initial 


squeezed coherent states in a quantum chaotic optical fiber. 
It is to be noted that a judicious choice of a quantum chaotic 
core cross-section would lead to the generation of highly en¬ 
tangled transverse modes. These results are also in accordance 
with our previous studies (0, IH] . 


V. CONCLUSIONS 

We have explored the dependence of the quantum entan¬ 
glement with the boundary deformation. We have seen that 
as the boundary become deformed into the quantum chaotic 
regime, the entanglement in the eigenmode increases. This 
demonstrates a simple approach to manipulate entanglement 
via mechanical means. We have also analyzed the propagation 
of a coherent state and a squeezed coherent state in a quan¬ 
tum chaotic Robnik optical fiber. We have found that an ini¬ 
tial squeezing can indeed enhance entanglement in quantum 
chaotic fibers. More importantly, our results have specifically 
proven the efficiency of manipulating entanglement in an op¬ 
tical fiber by modifying its cross-sectional geometry. 
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